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Abstract 



We study infrared QCD effects in radiative quarkonia decays. We examine the 
endpoint region z — > 1 of the photon spectrum. We point out a cancellation 
mechanism for the corrections in a"ln m (l — z), m < 2n, in the short-distance 
coefficient for the color-singlet Fock state in the quarkonium. The cancellation is 
due to the coherence of the color radiation, and applies even though logarithmic 
contributions are present in the jet distributions associated with the decay. We 
comment on the implications of our results for the modeling of hadronization 
in the endpoint region and for the role of color-octet states in the quarkonium. 
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1. Introduction 



Radiative decays of heavy quark-antiquark bound states (quarkonia) have been inves- 
tigated in the framework of perturbative QCD and have been used to measure the QCD 
coupling at energies of the order of the heavy quark mass m ]T|, ^ [| . The use of pertur- 
bative QCD in this context is based on the fact that, as long as the velocity v of the heavy 
quark is small, these decay processes involve two widely separated distance scales: the 
scale l/(mv 2 ) over which the quark and antiquark bind into the quarkonium and the scale 
1/m over which the quark-antiquark pair decays. By expanding about the nonrelativistic 
limit v — ► 0, one may treat the process as the product of a long-distance, nonperturbative 
factor, containing all of the bound-state dynamics, and a short-distance factor, describing 
the annihilation of the heavy quark pair and computable as a power series expansion in a s . 

However, near the exclusive boundary of the phase space, where the photon's energy 
approaches its kinematic limit, both the expansion to fixed order in v and the expansion 
to fixed order in a s become inadequate to represent correctly the physics of the decay. On 
one hand, classes of relativistic corrections that by power counting are higher order in v 
become enhanced by powers of l/a s || || [5J. On the other hand, potentially large terms 
in ln(l — z) (with 1 — z being the parameter that measures the distance from the endpoint) 
appear in the coefficients of the expansion in a s to all orders ||. In both cases, reliable 
results can only be obtained after resummation of the enhanced contributions. 

Infrared QCD effects are responsible for these behaviors. Soft color interactions account 
for the heavy quark-antiquark pair turning into a color-octet state, whose contribution, 
although subleading in v, becomes important near the endpoint. Sudakov logarithms in 
1 — z arise as remnants of the imperfect cancellation of the infrared divergences in the 
exclusive region, where particle production is suppressed. The need for resummations noted 
above is a symptom of the sensitivity of the process to contributions from infrared scales. A 
corresponding effect is seen in calculations that incorporate models for the hadronization of 
partons ||, ||, [K| , where nonperturbative contributions are found to be essential to describe 
the endpoint region of the photon energy spectrum in J /if) |TTJ and T || [L2| decays. 

In this paper we present a treatment of the infrared QCD radiation near the kine- 
matic boundary. We describe contributions from soft and collinear gluons with logarithmic 
accuracy, and study the color correlations between the jets that accompany the photon. 
This treatment takes account of the coherence properties of gluon emission. We determine 
the regions in the final-state phase space for the quarkonium decay in which destructive 
interference between parton emitters suppresses color radiation. 

We use these results to examine potentially large Sudakov logarithms in the photon 
spectrum. We find (see Sec. 4) that Sudakov logarithms cancel order by order in a s in 
the short distance coefficient for the color-singlet Fock state in the quarkonium. The 
cancellation mechanism applies to the leading logarithms as well as to any subleading 
logarithms. As a result, no Sudakov suppression factor arises near the endpoint. 
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The picture underlying this behavior can be understood in simple terms based on the 
branching formulas for the decay discussed in Sec. 2. In the boundary kinematics the pho- 
ton recoils against two almost-collinear gluon jets. The cancellation of Sudakov corrections 
reflects the fact that color is neutralized already at the level of this two-jet configuration. 
This situation may be contrasted with the situation one encounters in the decay of an elec- 
troweak gauge boson into jets. Here Sudakov corrections arise precisely from the presence 
of color charges in two-jet configurations. 

This picture also indicates that no cancellation should occur for the color-octet Fock 
state in the quarkonium. In this case, we expect the usual Sudakov suppression to take 
place near the endpoint of the photon spectrum. Then one of the consequences of the 
results of this paper concerns the ratio of the color-octet to the color-singlet contributions. 
This ratio will be smaller than expected from the power counting in v and a s obtained by 
truncating perturbation theory to fixed order || [13|, 14| , owing to the different high order 



behavior of the perturbation series in the two cases. 



A preliminary discussion of the results presented in this paper appeared in Ref. 15 



The paper is organized as follows. In Sec. 2 we treat the emission of soft and collinear 
gluons associated with the quarkonium decay. This leads us to a general branching formula. 
In Sec. 3 we evaluate the phase space in this formula. In Sec. 4 we compute the photon 
spectrum and observe the cancellation of the Sudakov logarithms. In Sec. 5 we discuss the 
implications of our results and give conclusions. 



2. Infrared radiation near the phase space boundary 

In this section we describe the QCD radiation that accompanies the decay of the quarko- 
nium near the endpoint of the photon spectrum. The key element in our treatment is the 
coherent emission of gluons and the angular ordering of this emission. We arrive at an 
expression for the photon spectrum that is valid to leading and next-to-leading accuracy 
in the soft and/or collinear logarithms. 

We consider the decay of a quarkonium state H of mass M into a photon plus anything: 

H^-f + X . (1) 

We concentrate on the case of the leading Fock state for a 3 Si quarkonium state, that is, 
a heavy quark- ant iquark pair in a color singlet state. 

This decay receives contributions at the leading order of perturbation theory both when 
the photon is directly coupled to the heavy quarks (direct term) and when the photon 
is produced by collinear emission from light quarks (fragmentation term) |Tj|. In the 
endpoint region the fragmentation contribution is suppressed relative to the direct term. 
The discussion of this paper will thus focus on the direct term. 
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2.1 Notations and tree-level decay 



The direct contribution at the leading order is given by the process 

H(P)^ 1 (k)+g(k 1 )+g(k 2 ) (2) 

evaluated at the tree level. We introduce dimensionless energy variables for the particles 
in the final state, as follows 

In the quarkonium rest frame z and are respectively the photon and gluon energies scaled 
by the heavy quark mass m (m ~ M/2 in the lowest nonrelativistic approximation). 



The normalized photon spectrum from the process (g) is given by [17] 

1 dT 1 r 1 r 1 

dxi / dx 2 p(xi,x 2} z) 5(z — 2 + xi + x 2 ) 0(xi + x 2 — 1) , (4) 



To 7T 2 — 9 Jo 

with 

P(^l»^.-g)= o T 2 + o r 2 + r 2 r 2 • ( 5 ) 

Over almost the entire range in z the result (f|) is well approximated by a spectrum rising 
linearly with z. At the endpoint, however, the first derivative of the spectrum (|4]) is singular. 
We will come back to this in Sec. 4. 



2.2 Behavior at z —>■ 1 

Beyond the tree level, corrections in ln(l — z) to the spectrum become possible |§. These 
come from the infrared region: near the endpoint the emission of gluons is suppressed and 
terms in ln(l — z) result from the imperfect cancellation of the infrared divergences between 
real and virtual graphs. The contribution from the soft and/or collinear gluons associated 
with the decay is at most double- logarithmic for each power of a s . By power counting the 
leading behavior of the photon spectrum as z — > 1 is of the type || 

- j- ~ const. +]Tc fc a k s ln 2fc (l - z) , z - 1 . (6) 

k=l 

In this paper we will see a cancellation mechanism for these corrections that works 
order by order in perturbation theory. A qualitative picture of the cancellation that we will 
discuss in the next sections can be given at one loop as follows. Consider the emission of 
a soft gluon with momentum q from the tree- level graphs for the decay @. In the leading 
infrared approximation, the amplitude A for this process factorizes (see for instance |T%fl) 
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into a nonabelian, classical current j(q), describing the soft-gluon emission, times the tree 
level amplitude Aq- 



A b (P, h, k 2 - q) ~ g s e»(q)jl(q)MP, h, h) , (7) 
where e is the gluon polarization vector, b is the gluon color index, and 

Here the sum runs over the final-state colored particles k\ and k 2 , and Ti are the color 
matrices associated with the coupling of the soft gluon to each of these particles. There 
is no contribution from the coupling to quarks: the coupling to the virtual quark lines is 
subleading in the infrared limit because these lines are off shell by order M; the coupling 
to the quark and antiquark in the quarkonium does not contribute either, although these 
quarks are nearly on shell, because in the lowest nonrelativistic approximation they are 
almost collinear and their total color charge is zero. 



To the approximation at which we are working, the phase space also factorizes |18 
that one gets a factorized answer for the spectrum, schematically of the form 

dT r - dT 



a 



Here d<&' is the Lorentz-invariant phase space for the soft gluon emission. The standard 
power counting in terms of the soft gluon energy uj gives 

dV-ududn, i 2 ~— 2 fl kl ' k l , , (io) 

or j{z; angles) 

where dfl is the angular phase space, and / is a function of z and the angles as uj — > 0. Up 
to the first order in uj the correlation of the gluon momenta is 



k x -k 2 ~ M 2 



z) + J^9{z; angles) 



(11) 



That is, as z — > 1 the photon recoils against two almost-collinear hard gluons. In this 
configuration, the logarithmic integration du/uin Eq. (|9|) is canceled. The next few sections 
are devoted to showing that this is indeed the mechanism that dominates the radiation 
accompanying the decay, and that it extends to all orders in a s . 

The back-to-back kinematics suggests an analogy with the two-jet region in the anni- 
hilation of e + e~ into electroweak gauge bosons. Since the two jets are now in a singlet, 
though, the pattern of color correlations will be different. Jet event shapes in e + e~ annihi- 
lation have been studied extensively by using methods based on branching graphs. We will 
see in the remainder of this section that the photon spectrum in decays of quarkonia can 
be treated in a similar manner. The subtlest point will be the treatment of the branching 
phase space. This will be addressed in Sec. 3. 
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2. 3 Branching graphs 



Methods to calculate multi-parton matrix elements in the leading soft and collinear 
orders have been known for a long time. See Ref. jy| for a review. The essential observation 
in these methods is that one can replace the calculation of higher-loop Feynman graphs by 
the calculation of angular-ordered "branching" graphs, that is, tree level graphs in which 
the angular phase space is subject to certain ordering constraints at each branching. The 
method has also been extended to include next-to-leading logarithms: see for instance 
Ref. |19|] . In this and the next subsection we apply this approach to the quarkonium decay. 



We write the energy spectrum dT/dz in the inclusive decay (JJ) in terms of the exclusive 
decay widths dT^ xc ^ for producing 7 + n final state partons qi, . . . , q n : 



(12) 



Here the 5 function expresses the photon energy fraction z in terms of the quarkonium and 



parton momenta. The exclusive decay widths are generated by the branching process ||T8 

Eg. 



The first step of the branching corresponds to the lowest order final state, Eq. (Q), 
in which the photon is accompanied by two gluons with invariant mass (ki + k 2 ) 2 . This 
invariant mass provides the basic mass scale at the hard end of the branching. The sub- 
sequent steps correspond to consecutive parton splittings (Fig. la). There are well-defined 
probability weights associated with each splitting vertex, and form factors associated with 
each line connecting two vertices [|T^, |TDJ . The basic point about the branching process is 
that the phase space for the splittings is restricted to the angular ordered region in which 
the branching angles decrease as we go from the parent parton toward the final state. If we 
denote, as in Fig. la, by i/i the energy transfer at the vertex i, and by r^ the transverse 
momentum flowing between vertices i and i + 1, this region is defined (for small yi) by 



(kt + k 2 



vl vl 



> 



(13) 



The reason for the angular ordering is a coherence effect: outside the ordered region parton 
emitters interfere destructively pO , 21], 



]. In this approach, the exclusive decay width 
^p(exci) j g obtained by taking all angular-ordered branching graphs with n final state par- 
tons, and summing the corresponding probability weights. 



Following |19| , we integrate the exclusive decay widths over the phase space of the final 
state partons keeping k\ and k 2 fixed. We rewrite the n-parton phase space as (Fig. lb) 



- r/V 



d%5 4 [h-J^Qj 



x d 4 k 2 5 A k 2 




1 (2vr) : 



n d A n- 

n tSm^ 



(14) 
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Figure 1: (a) Branching process; (b) final states in the quarkonium decay. 



The integration of the branching matrix element over the final state phase space at fixed 
k 2 (i = 1,2) produces the jet mass distribution J g (p 2 ,k 2 ) |L9). This is defined as the 



probability to generate a jet with invariant mass k 2 from a parent gluon produced at the 
mass scale p 2 . We have 



J 9 ({ki + k 2 )\kl 



1 d A n- 

n tSrMfi 



f=\ (2vr) 3 
(branching matrix element) 



(2tt) 4 5 4 h 




(15) 



and an analogous expression for J g ((ki + k 2 ) 2 , k\). Note that the jet mass distributions 
depend on the mass scale ik\ + k2) 2 . As implied by the coherence relation (13), this is the 
mass scale at which the initial gluons in the branching are produced. 



2.4 Jet mass distributions and inclusive spectrum 



To lowest order, there is only one final state parton in the definition (y|) of the jet mass 
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distribution, and J g is just given by a delta function: 

J g [p\ k 2 ) = 5{k 2 ) + 0{a s ) 



(16) 



In general, the analysis of the jet branching process allows one to obtain an evolution 



equation [19| for the jet mass distribution. Schematically, this has the form 

J(p 2 , k 2 ) = 5(k 2 ) + J a s (p' 2 )IC(p' 2 , k 2 ) ® J(p' 2 , k 2 ) , (17) 
where K, is a kernel calculable as a power series expansion in a s . Expressions for the kernel 



/C are known to the leading [£L8|, [2(], ^] and next-to-leading order [19], |23|]. Solving the 
evolution equation to the leading logarithms gives 



J dk 2 J g (p 2 , k 2 ) &(K 2 - k 2 ) = exp In (j^j f (f3 a s In 

where (3$ = (HCa — 2iVy)/(127r) and / is the function 

C A 1 



K 2 



2nP x 



[(1 - 2x) ln(l - 2x) - 2 (1 - x) ln(l - x)\ 



18) 



(19) 



It is also useful to introduce the double logarithmic approximation. This is defined by 
expanding the exponent in the right hand side of Eq. flT8| ) to the first order in a s . In this 
approximation, using f(x) = —Cax/(2tiP ) + 0(x 2 ), we get 



J dk 2 J g (p 2 , k 2 ) Q(K 2 - k 2 ) w exp 



(20) 



In order to express the inclusive spectrum in terms of the jet mass distributions, we 
now use the definitions (14]), (15) in Eq. (|T^) . By taking account of the photon phase space 
factor, the inclusive spectrum takes the form 



1 dY 

r dz 



d A k d 4 k\ d 4 k 2 , n , , , s . / 2P ■ k\ . ,, 9 . 

:(2tt) 4 5 4 (P — k — ki — k 2 )5 z -— 5 + (k 2 ) 



(2tt) 3 (2vr) 3 (2vr) 3 



M 2 



x Mq(P, ki, k 2 ) J g Uki + k 2 ) , k{) JJih + h) 2 ^, 



(21) 



Here A4 is the amplitude for the first step in the branching, with final state ^ygg. Owing 
to the constraint set by coherence, we can use the approximation (£4 + k 2 ) 2 ^> k\,k\ and 
evaluate A4 by setting k\ = k\ = 0. The explicit expression is readily computed and reads 



Mc 



8vr 3 M 2 



7T 
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±{ki ■ k 2 ) 2 



+ 



(p - hY (p - k 2 ) 



(p ■ h) 2 (p ■ k 2 y (p-k) 2 \(p-k 2 ) 2 (p-h) 2 



(22) 



We are now going to evaluate the branching formula (pT|) explicitly. To this end, we 
need to analyze the phase space of this formula in detail. This is the object of the next 
section. 
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3. Structure of the phase space for small jet masses 



We are interested in the angular-ordered, coherent region 

k\ , k\ < (ki + k 2 f < M 2 . (23) 

Consider the phase space element in Eq. fl2~T|): 

rf$ = d 4 k d% d A k 2 6 4 {P -k-k 1 -k 2 )diz- 5+{k 2 ) . (24) 

It is convenient to express this phase space in terms of the energy fractions x\, x 2 defined 
in Eq. (£J), the jet masses kf, k 2 , and 

sn = (ki + k 2 ) 2 . (25) 

We can use the momentum conserving S in Eq. (^) to integrate over k. From the 
positivity of the jet energies and invariant masses we have 

0<xi<l, 0<x 2 <l, xi + x 2 >l , (26) 

< k\ < x\M 2 /A , < k\ < x 2 2 M 2 /4 . (27) 
The relation (^5|) between s\ 2 and the jet momenta implies 

k\ + k\ + l -x x x 2 M 2 - 2|ki||k 2 | < sia < k\ + k\ + \x x x 2 M 2 + 2|ki||k 2 | (28) 



with |ki| = ^x\M 2 /A - kl |k 2 | = ^x|M 2 /4 - ki 
The phase space element can then be rewritten as 

= _ M dxi dx 2 dk\ dk\ ds l2 5(z - 2 + x x + x 2 ) S(M 2 (1 - z) - s 12 ) , (29) 



with the constraints (p6D,(|27|). 



Eq. (^) selects in general a subspace of a complicated form. But we need to consider 
it in the logarithmic region (^). We may use the 6 functions in Eq. (|29|) to integrate over 
Si 2 and one of the energy fractions, say, x 2 . Then Eq. (|28f) gives 

x± < X\ < x\ , (30) 

with 

xf = (1 - z/2) (l-z + kj/M 2 - k 2 2 /M 2 )/(l - z) 

± {z/2) \j(k 2 - k 2 ) 2 /M 4 + (1 - zf - 2 (1 - z) {k\ + kl) /M 2 /(l - z) . (31) 
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By approximating the constraints (|30|) for kf, k\ <S M 2 (l — z) <C M 2 , we get 

■ Tl - M 2 (l-,) ' (32) 

and 

^-'"M^) • ( 33 ) 



that is, since 1 — xi ~ X2 for 2 — ► 1, 



^M^Ti • (34) 



Eqs. fl32|) and fl3~4|) tell us that in the coherent region the phase space available for the 
evolution of each of the jets is bounded by the recoiling jet. For z — > 1 this bound is tighter 
than the bound from the fragmentation of the jet itself, Eq. (0). In the next section we 
will see that the recoil constraint conspires with the angular-ordered form of the jet mass 
distribution to destroy any logarithmic hierarchy in the inclusive photon spectrum. 



4. The photon spectrum near the endpoint 

We now put together the results of Sec. 3 for the phase space with the structure of the 
spectrum described in Sec. 2 based on the coherent branching. 

By using Eqs. (j22|) , (p9"|) , (|32]) and fl34| ) in Eq. (plj), we can rewrite the photon energy 
spectrum as follows: 

1 dT lr 1 r 1 

— ~r — —5 — - / dxi / dx 2 p(xi,x 2 , z) 8(z — 2 + xi + x 2 ) Q(xi + x 2 — 1) 
1 dz 7r — 9 jo jo 

x f ikl J, (M^ - ,), *») 9(^(1 - ,) - k\) - kl) 

x / dk 2 2 J g (M 2 (l - z), fcfl 9(M 2 x 2 (l -z)- kj) Q(M 2 xl/A - k\) , (35) 

<J 

with p(xi,x 2 ,z) given in Eq. (||). The main structural difference with respect to the case 
of jet event shapes in e + e~ annihilation is that in Eq. (]35f) there is a two-dimensional 
integration over a distribution p in the energy fractions x\, x 2 . This comes from the fact 
that in quarkonia decays the parton branching is probed by a non-pointlike source. Eq. (j35|) 
allows us to discuss the logarithmic behaviors in the endpoint region. 

Observe first that, by substituting into Eq. ( |35"D the zeroth-order expression (0) for 
the jet mass distribution J g , we recover the lowest order result (f|). So at the lowest order 
of perturbation theory Eq. (^) gives the correct answer for any z. At higher orders of 
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perturbation theory, Eq. (|35|) gives an approximation valid in the region of large z with 
leading-logarithm and next-to-leading-logarithm accuracy, provided the jet mass distribu- 
tions J g are evaluated with corresponding accuracy. Once expanded to the next to lowest 
order in a s , Eq. (^) can be matched with the NLO perturbation theory result [2l|, and 
could thus be used to obtain improved predictions, valid over a wider range of z. 



Eq. ( p5| ) enables us to see that, although Sudakov logarithms are present in the jet 
distributions associated with the decay, the corrections in ln(l — z) to the photon energy 
spectrum cancel order by order in a s . The mechanism for the cancellation can be seen in 
its simplest form by using the double- logarithmic approximation (|20D for J g . By expanding 
the right hand side of Eq. ( GDp in powers of a s and substituting this into Eq. (^), we note 
that the higher order corrections to dT /dz involve integrals of the form 

That is, logarithmic contributions in x\ arise, which are important at the kinematic limit 
x\ — > 0, but these never give rise to logarithms of (1 — z) in the photon spectrum. Terms 
in ln(l — z) cancel because of coherence, 1.6. ? clS db result of the constraint on the jet 
mass and of the angular ordering (|13D in the branching. 



From Eq. fl35|) we can obtain an explicit expression for dT / dz by keeping track of all the 
leading logarithms in the jet distributions. This is accomplished through Eq. (|i~8f) . Using 
this formula in Eq. fl35|) we get 



1 dT 



/ dx\ \ dx2 p(xi,X2, z) 5(z — 2 + x\ + X2) 0(xi + X2 — 1) 
Jo Jo 



r dz 7T 2 - 9 

x exp KIM) / In(lM)) + \n(l/x 2 ) f ((3 a s ln(l/x 2 ))] , (37) 

where f(x) is given in Eq. (^). It is straightforward to check from Eq. (|37|) , by expanding 
the integrand in powers of a s , that no logarithms of (1 — z) appear in the perturbative 
expansion of dT /dz. 

The above results indicate that the photon spectrum from the decay of the color-singlet 
Fock state in the quarkonium is not Sudakov suppressed by higher orders of perturbation 
theory in the endpoint region. Higher perturbative corrections give rise to a constant shift 
compared to the leading order answer (|J). 

The first derivative of the spectrum, on the other hand, does get logarithmic corrections 
in (1 — z). This can be realized just from Eq. (|36|) by taking the derivative with respect 
to z. A singularity in the derivative of the spectrum at z — 1 was indeed noted already 
in lowest order (see comment below Eq.(||)). The leading logarithmic behavior at one loop 
will be of the type 



d (I dT s 
dz \ T dz 



a In (l-z) + const. + 0(1 - z) 1 } (38 



+ oc. 



at In 3 (1 - z) + subleading logs + 0(1 - z) 1 ] + 0(a]) 
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It is of much interest to calculate these corrections. Likely, the z = 1 divergence that 
appears in the derivative of the spectrum at any fixed order of perturbation theory will 
be smoothed out by the all-order summation of the logarithms. Note however that this 
calculation will not necessarily involve only the configurations in which the photon and the 
gluon jets are at large relative angle, and the jets evolve according to the angular-ordered 
branching, but also configurations in which color is emitted at angles comparable to that 
of the photon. 

To conclude this section, we remark that near the boundary the decay becomes sensitive 
to the hadronization process. For the inclusive spectrum, these effects can be factorized 
in nonperturbative shape functions |2J| |27|. A very important question concerns the 
modeling of these functions. The Monte Carlo calculation || provides one such model, 
based on parton showering and the assumption of independent fragmentation. Another 
model is that of Refs. |4], 10| , in which nonperturbative corrections are parameterized in 
terms of an effective gluon mass m g | 2q| . We observe that results on the infrared behavior 
such as those presented in this paper can be used to investigate models of hadronization. In 
particular, if the resummed formulas for the photon spectrum discussed above are combined 

they 



with models for the behavior of the running coupling at low energy scales [^6|, |29| 
provide an ansatz for the power corrections in the shape functions [|30|]. These power 
corrections are the analogues of the contributions in m g /M of [TO, 
to be necessary to describe the data in the endpoint region. 



23, which are found 



5. Conclusions 

The photon spectrum in quarkonia decays is a critical issue in QCD phenomenology. 
Although, as was realized early on ||31|| , for large enough masses the decay is dominated 



by short distances, the observed behavior [y, |TT|, |TI| of the spectrum at large z is not 
reproduced by fixed-order perturbation theory. In this respect the first implication of the 
results of this paper is negative: perturbative resummations do not fix the problem. We 
have seen in the previous sections that resummation does not produce a large- z suppression 
of the spectrum, but a constant shift. 

A second, perhaps more general implication is that color coherence effects in the quarko- 
nium decay are important. As we have seen, they change dramatically the large- z behavior 
of the perturbation series compared to what one would conclude from the simple infrared 
power counting ||, leading to the cancellation of the Sudakov terms. 

Note that the comparison of theory with experiment (see, e.g., 0, [3|]) has so far involved 
the use of models for the parton cascade and hadronization |5] that do not include coherence. 
In the case of |J the emission of color associated with the evolution of the jets takes place 
within a cone whose typical size is 9 ~ 1. But the analysis of this paper indicates that in 
fact, since the coherence scale in the jet mass distributions is not the quarkonium mass M 2 
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but rather {k\ + k 2 ) 2 , destructive interference occurs outside a cone with size 9 ~ \/l — z. 
It will be valuable to make Monte Carlo models for quarkonium decays that take this into 
account. 

In the endpoint region relativistic corrections may become important |5|, 13, 53|. De- 
tailed estimates of color octet contributions have recently appeared || [jj. In this paper we 
have observed the absence of Sudakov suppression for the Fock state of lowest order in the 
nonrelativistic expansion, i.e., a color singlet quark- ant iquark pair. The mechanism that 
produces this result is based on the color correlations of two gluon jets, and therefore is 
not at work in the case of the color-octet state. This means that the power counting in the 
heavy quark velocity and in a s at fixed order is modified by the high order behavior in very 
different ways for the color-singlet and color-octet channels. The relative size of the color 
octet with respect to the color singlet in the endpoint region will be smaller than indicated 
by the power counting at fixed order. 

Quarkonia decays are used to measure the QCD coupling and contribute significantly to 



the world average value of a s |$3|, |4j . In fact Ref . |34j quotes a very precise determination 
of a s using this method. As we have seen, however, effects from the endpoint region in 
radiative decays may be quite dramatic (cancellations at large z, lack of angular ordering 
in current Monte Carlo models, need for an all-order power counting for relativistic correc- 
tions), and have yet to be fully understood. The present uncertainty on the extraction of 
a s from these decays is therefore likely to be larger than previously estimated. 

There has recently been progress in next-to-leading-order calculations for the photon 
spectrum. NLO corrections have been computed for color octet contributions || and for 
color singlet, direct contributions ||24|| . The missing piece are the corrections to the frag- 



mentation terms (JlfJ, to be used with the next-to-leading fragmentation functions of the 
photon [35]. A possible choice for future determinations of a s will be to use only data 
sufficiently away from the endpoint for NLO perturbation theory to be a valid approxi- 
mation 10, This raises the question, though, of where the safe region begins and how 
much of the data one is left with. Alternatively, note that resummation formulas of the 
type presented in this paper can be matched and combined with NLO results. This will 
give improved predictions, which could likely be used in a wider range of photon energies 
towards the peak region. 

Hadronization physics, showing up as power-suppressed corrections to the spectrum, 
becomes a dominant effect near the endpoint. Power corrections are indeed found to be 



essential to describe the data in this region |4J, [T(J. It becomes important to investigate 
models for the nonperturbative shape functions |26|, parameterizing these effects. 

The structure of soft color emission discussed in this paper, once combined with an ansatz 
for the infrared behavior of the coupling |30| , may serve to study these models. 
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